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Abstract 

We study the Schrodinger-Newton system of equations with the addition of gravitational field 
energy sourcing - such additional nonlinearity is to be expected from a theory of gravity (like 
general relativity), and its appearance in this simplified scalar setting (one of Einstein’s precursors 
to general relativity) leads to significant changes in the spectrum of the self-gravitating theory. 
Using an iterative technique, we compare the mass dependence of the ground state energies of 
both Schrodinger-Newton and the new, self-sourced system and find that they are dramatically 
different. The Bohr method approach from old quantization provides a qualitative description of 
the difference, which comes from the additional nonlinearity introduced in the self-sourced case. 
In addition to comparison of ground state energies, we calculate the transition energy between the 
ground state and first excited state to compare emission frequencies between Schrodinger-Newton 
and the self-coupled scalar case. 
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I. INTRODUCTION 


The Schrodinger-Newton (SN) system of equations has been studied in various contexts 
since its introduction in pQ as a model for self-gravitating (quantum) particles. The single¬ 
particle formulation of the problem treats the central body’s quantum mechanical distribu¬ 
tion as a mass density sourcing the gravitational field (whose influence affects the central 
body): 

,,9\l/(r,f) X NT/ X 

ih ---=-V^\I/(r,f) -I- m$(r,f) T(r,f) 

ot 2 m 

V^<h(r,t) = 4 7rGmT*(r,f) T(r,t) 

where T(r,f) is the wave function associated with the particle of mass m moving under the 
influence of the potential energy $(r,f). 

Regardless of its motivation (and there are a variety of motivations for studying this set 
- some of these are described in [2[ |T9] and references therein - we side-step that discussion 
in the current work, focusing instead on a comparative study of the solutions to ([^ and our 
proposed modihcation), much is known about the solutions in both the time-independent [3], 
m and time-dependent form [5HZ] (see also the dissertations 0191 for additional background 
and review). 

We propose to augment the gravitational sourcing in 0 to include the energy density 
present in the gravitational field itself. The motivation for doing this comes originally from 
special relativity and the universal coupling of gravity - if mass density acts as a source, 
then so can energy density, and the energy density of the gravitational field is available 
as a source prior to any additional external sources (like the energy density associated 
with electromagnetic fields, for example). The governing equation for $, the gravitational 
potential energy, is then [10] : 

V^4(r) = ^ p(r) 4>(r) + V4(r) • V4(r), (2) 

2<h(r) 

This is the held equation for a scalar <f)(r) that is sourced by a distribution of mass given 
by the mass density p{r) and its own energy density (represented by the second term). This 
held equation was developed originally by Einstein [lO], and a (special) relativistic version 
appeared in m - more recently, the self-coupled case was re-derived in [T2] . 

If one were to similarly couple a second-rank held theory to itself by making the held’s 
stress tensor a source - motivated by the universal nature of gravity - one ends up with 
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general relativity (GR) [13] (and for a review of that process [HI) together with its geometric 
interpretation. The sourcing in (|^ is a scalar version of that self-coupling (and was worked 
out in that context in [H]). Our fundamental question will be how this relatively simple 
nonlinearity, showing up in the gravitational held equation [? ] changes the spectrum of the 
resulting Schrodinger-coupled system. 

We’ll start with a review of the self-coupling that leads to the (static) held theory (|^, 
then develop the spherically symmetric vacuum solution of ([^, and use that to estimate 
the ground state energy as a function of particle mass. Then we will introduce a numerical 
method to calculate the ground state energy of both SN (to test correct numerical behavior) 
and the new self-coupled gravitational held, and compare the ground state energy to the 
predicted form. For both the ground state and the hrst excited state, the Bohr model 
estimates provide qualitatively correct predictions. By calculating the energy diherence 
between the hrst excited state and ground state of the self-coupled gravity, we can compare 
the transition energy emission of this theory with the one predicted by SN. Finally, we 
calculate the total energy of N self-gravitating bosons using the self-coupled form, and 
compare that with the total energy as calculated using SN in [1]. 

II. MOTIVATION 

The static gravitational held theory given by (|^ was hrst proposed by Einstein en route 
to general relativity. We will generate the held equation here (a recent discussion of this 
derivation is in pEI) - the issue that motivates its development is the lack of self-coupling of 
the original held equation of Newtonian gravity: = AnGp describes a held $ sourced 

by a mass density p, but this Poisson form lacks the energy density source that comes from 
<h itself. 

Let’s start by calculating the energy density associated with a held $ that comes from 
Newtonian gravity (so that it satishes the Poisson equation: = An G p). Proceeding as 

usual, the work done in building a distribution of mass, p, is given by: 

f p<S>dT, (3) 

all space 

then using the Poisson equation to eliminate p and integrating by parts gives 
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and we would call the integrand the energy density of the held: 

= (5) 

O TT (jr 

According to special relativity, we should use an associated effective mass density 

as a source in the held equation for $. That is the nonlinearity promised by a theory of 

gravity - the held’s self-energy must act as a source. So we are tempted to start with: 

V^d) = 47rG'(p + p$) =47rGp-^V<h-V<h, (6) 

2 & 

and this form was considered in HU. The problem is that the new held equation ([^ leads to 
an energy density (obtained as above, but with the new held equation instead of the Poisson 
equation) that is not (|^, and so the self-coupling is not consistent. 

We’ll now generate the correct self-coupled theory using the machinery of held theory, 
where the Hamiltonian density is the energy density (for a static, free held). Suppose we 
start with a free held theory that has Hamiltonian "H = V<h ■ V4>, i.e. we augment the 
usual scalar Hamiltonian with a function of <h, /(*h), that we will hx by demanding that it 
lead to a source term that looks like = "H/c^. The held equation for this "H is 

ZM v^d) + V4> • V<h = 0, (7) 

47rG 8nG 


or 


V^<h = 


/'(«>) 


2^(4)V4.V4. (8) 

It is the right-hand side of this equation that we would like to set equal to dvrG (Tf/c^) = 
V<h ■ V<h, giving us the self-consistent energy density source, and when we do that, we 
get an ODE for /(<h): 

/($) /'($) 


fW = 


( 9 ) 


2c^ 2 /(®) 

2 _ _ 

Our hnal free Hamiltonian is: "H = V4> ■ V<h, and this is the energy density of the 

held in the theory, given also in [10]. Using this solution for /($) gives the held equation: 


V^<h = — V<h ■ V<h, 
2 <l> 


( 10 ) 


in vacuum, and we recover (|^ when we introduce massive sources (that would put a factor 
of p on the right-hand side of ([^). When computing solutions to this self-consistent, self- 
coupled form of gravity, it is worth noting that the held equation ([^ can be written linearly 
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(11) 


in a/$, where it reads (again, from [10]): 

In order to recover Newtonian gravity, we must take $ = + <hAr, then inserting this 

into 0 (or (|^) and collecting in powers of c gives back V^<hAr = dvrCp to zeroth order, 
a requirement of the weak-held limit. 


III. POINT SOURCE SOLUTION 


We’ll start by looking at solutions to the held equation: 

+ — Vd)-V$, (12) 

2 d) 

in regions where p = 0 (so we are in vacuum) and $(r) = $(r). Such a solution would be 
appropriate for a spherically symmetric source of mass m localized near the origin. Under 
our assumptions, the held equation reduces to 


(r $(r))^^ = 


2$(r) 


(4'(r))'- 


(13) 


with primes denoting r-derivatives. The solution to this equation comes with two integration 
constants, a and (3: 


<I)(r) = (3 


. 9 4a 1 
4a —h — 


(14) 


If we ask that <I)(r) look like a Newtonian point source (of mass m) as r approaches spatial 
inhnity, then we can hx a 

Cm R m^ 

(15) 


, Gm (3 G‘^m? 

$(r) =-+ Ut + 


r 4/3 


To recover the Poisson form of Newtonian gravity for 4) small (compared to c^), we must 
have <I>(cx)) = as a boundary condition, and that sets the constant (3. Our hnal spherically 
symmetric vacuum solution looks like 

Gm G'^ 


(|)(r) =-h 

r 

The energy density of this held is, from 73 = 

7i = 


4 


+ cl 


(16) 


8 TT G $ 

Gm^ 

8 TT ’ 


V<I>, 


(17) 
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which is everywhere positive. 

We are interested in the vacnnm solntion becanse it is the one that is relevant to the Bohr 
approach to estimating qnantnm mechanical energies. Bnt it is easy to compnte (especially 
from (pBj)) the solntion in cases other than vacnnm - for example, if we had a sphere of 
radins R with constant density po inside it, then the interior solntion to (|^ is just: 

A sinh (r/ro) \ ^ c 


$i(r) = 


ro = 


(18) 


r J y/2 TT G Po 

where A is a constant that we would use to match up to an exterior solution (at r > R, for 
example) and we have chosen the solution that is finite at the origin. 


IV. SCALING OF THE GROUND STATE 


For the SN system, with constants G, h, and m, there is only one way to make an energy 

G2 




(19) 


SO we expect, up to constants out front, that the energy spectrum of ([^ scales like nR (as 

indeed it does mm)- 

With the introduction of c appearing in (2), we can form the Planck mass: Mp = 
and this means that any power of m could appear in the energy spectrum of the self-coupled 
system (where we expect E ^ with q + s = 5, but otherwise unconstrained). We’d 

like a way to estimate relevant combinations of Mp and m that might appear in our new 
spectrum. To that end, we will use the Bohr model (originally for hydrogen, of course, but 


applied in this gravitational setting) with potential given by (16) -- the idea is that if the 
ground state is localized close to the origin, then far away, the potential associated with that 


ground state should go roughly like (16), and so the spectrum of the spherically symmetric 


vacuum solution could provide some relevant approximate information. 

To start, we’ll apply the Bohr method to <l)(r) = — -f- c^, just the Newtonian point 

particle potential (with an offset at spatial infinity so as to match (Q). According to 
the rules of old quantization (for circular orbits, updated to elliptical orbits in the Wilson- 
Sommerfeld formulation), we start with the total energy for a particle moving in a circle of 
radius r (so that v = -s/Gm/r)-. 


E = me - 


1 Gm^ 


( 20 ) 
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Next, we assume angular momentum is quantized: L = mvr = nh for integer n, giving us a 


value for the radius: r = 


71^ h? 

G 


then using this radius in E we get a discrete set of energies: 


E„ = m c — - 


1 


2 77,2 h? 


( 21 ) 


The ground state corresponds to n = 1: 


z:. 2 1^2 5 2 1 5 

Ei = mc =mc 


( 22 ) 


which is what we expect, namely a linear (in m) offset (associated with the shift at spatial 
inhnity) and scaling. 


Performing the same procedure for the “point potential” in (16) gives 


2mc^ 

E„, = ^ 


+ 2 Mp ’ 


(23) 


with ground state energy: 


2mc^ 

Ei= ^ 


m^ + 2Mf 


(24) 


where again, we only care about the mass scaling here - our estimate cannot predict constant 
offsets and/or overall constants out front. We’ll come back to those later on. Note that the 


expression (24) reduces to (22) in the m Mp limit, as it should: 


2mc^ I ^ (f m 


m c 


(25) 


for m small. 


V. NUMERICAL APPROACH 

We’ll start by specializing to spherically symmetric solutions, then we’ll render the equa¬ 
tions of interest dimensionless, and put them in time-independent form to dehne the eigen¬ 
value problem of interest. From there, we’ll introduce the iterative hnite difference approach 
that can be applied to either SN or self-coupled gravity to hnd the ground state energies. 
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A. Dimensionless Form 


Since we are interested in the ground state energies, we will focus on spherically symmetric 
solutions to both SN and the modihed system. For SN, we have 


-(r'F(r, t)) + m <h(r, t) \h(r, t) = ? h 


2m r 


dt 

X 

- (r $(r,t))''= 47rGm'l/(r,t)* 


(26) 


and we can simplify further by setting P{r,t) = r4/(r, f), then the above becomes 


P(r, t)" + m <h(r, t) P{r, t) = ih 


2m 


(r$(r, t)) = 


dt 

„ An Gm 


(27) 


r 


P(r, t)* P{r, t). 


The wave function is normalized to 1, so that our P{r,t) has: 

poo 

An / ||P(r,t)||^ dr = 1. 

Jo 

The modihed theory becomes, under the same assumptions and substitutions, 
P(r, t)" + m <h(r, t) P{r, t) = ih 


(28) 


2 m 


(r<F(r, t)) = 


dt 

„ An Gm 


r 


P{r, t)* P{r, t) <h(r, t) + 


2 <h(r, t) 


(<h(r,t) 


/n2 


(29) 


We can render the equations dimensionless by introducing r = tq R, t = toT for dimen¬ 
sionless R and T, and taking <F = $o <5, P = Pq P, where $0 is a speed^ and Pq has dimension 
l/-\/length. Finally, let m = mofh for Planck mass mo = and set P{r,t) = e“*'®^P(r) 
for dimensionless energy E = E/Eq. Then our pairs take the form of an eigenvalue problem 

^ ^ ^ +m<lP = EP 

(30) 

P* P. 

R 

and 


m (9P2 


(R^) = — p* p 

dR2 


1 d^P 


m 


dR^ 


m^ P = E P 


(i? $) = — p* p $ H-- 7— 

^ ' R 2^\dR 


with 


h 


ro = 


V2 


tf) — 


dR^ 


h 


Pn = 


mo c 


mo P 


\/AnmQ G 


$0 = 


(31) 


Ea = nil, c'n 


(32) 



The normalization of the wave function (28) now reads: 


P*PdR = 


'/2Gm 


4 7rP2 


0^0 


he 


= V2. 


(33) 


B. Method 


The numerical method is the same in both cases - we discretize in R by taking Rj = j AR, 
where j = 0 is a boundary point (at the origin - for T hnite at the origin, we must have P = 0 
there), and we take Pj+i = Poo to be a numerical approximation to inhnity, where we again 
require P = 0. Let Pj = P{Rj) (and = |)(Pj)), then we can discretize Schrodinger’s 
equation using hnite differences: 

+ m^jP, = EPj (34) 


1 

rh 


Pj+i - 2 Pj + Pj-i 

AP2 


for j = 1... J. We can dehne the vector P G to have entries that are precisely the 


unknown Pj values (and similarly for the vector $), and then (34) can be written as a 
matrix eigenvalue problem in the usual way: 


D($)P = PP 


(35) 


with 


/ 2 
m A_R 2 


+ m 01 


1 

m A_R 2 


0 



D($) = 


1 

fh Ai?2 


2 

fh A_R2 


+ m 02 


1 

fhAR'^ 


0 


(36) 


V ° ■■■ ■■■ 'J 

As a matrix eigenvalue problem, it is relatively easy to construct D and then hnd the 
eigenvector associated with the smallest eigenvalue - that eigenvector is an approximation 
to the ground state. Once we have P in hand, we can construct the entries of $ in either 
the SN or augmented case using Verlet. We can hnd the values for by working backwards 
from j = J to 1 using the recursion 


$,_i = 


Rj-i 


2 Rj - P 


j+i <hj+i 


m 


+ ^\\pjr^R" 

ltd 


(37) 


where we take = 1 — for j = J and J + 1 - since Pj+i = Poo, we want the 
potential to approximate its value out at spatial inhnity, and this is the dimensionless form 
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of = c^ — For the gravitational field equation in (|31|), the analogous Verlet recursion 
looks like 


$ ■ 1 = 
^. 7-1 


Rj-i 


^3 -‘J -V+1 ^ {^3+1 ~ *^i) 


2 - Ri.-i <h,. 


(38) 


with the same boundary conditions as above. 

So, if we had the entires of P, we could construct the entries of $ in either case, but in 


order to get P, we need $ - the matrix D in (36) depends on the values of the potential $ 


We can use an iterative approach (similar to the one in [18]) to get around the problem. Let 


Pj be the iteration for Pj, then we can construct using (37), and update by finding 
the smallest eigenvalue/vector of D(<f)^), so 

P^^^ = smallest eigenvector of D($*^) normalized so that 'YTj=i ^R = V^- (39) 

From this, we can construct and iterate until: 


)fc+i 


P II < e, 


(40) 


for e, some user-specihed tolerance. 


VI. ENERGIES 


A. Schrddinger Newton 


The ground state energies for the SN system were calculated for m = .7 to m = 2 in steps 
of .02. In Figure we show P and $ for the masses rh = .7 (top) and m = 2 (bottom) - 
we chose these mass limits because at our value of Poo = 200 and AP = 21^^ masses less 
than rh = .7 begin to violate our assumption that P(oo) = 0, and masses larger than m = 2 
get localized to only a few grid points near the origin. We can expand the mass range by 
changing Poo, allowing us to probe smaller masses, and by decreasing AP, allowing us to 
move up in mass. 

The energies themselves are shown in Figure There, the dots are the numerically- 


determined values - we £t the dots to a curve of the form Am + B m^, motivated by (22), 
using a nonlinear Levenberg-Marquardt fit of the data. The curve is: 




C{m) = 1.00123m - 0.163181m 


(41) 
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FIG. 1. The numerical solutions for P and <I> for m = .7 (top) and m 
system. Here, i?oo = 200 with Ai? = 200/5001, and we use e = 10“' 


iteration (see (40)). 


= 2 (bottom) for 
as the tolerance 


the SN 
for the 


El 



FIG. 2. Ground state energy as a function of mass for the SN pair - the points are numerically 
determined using the numerical method described in the previous section, and the curve is the best 


fit curve, C{m), from (41) 


The best fit curve has two interesting features - first, it correctly identifies the linear mc^ 


offset from (22), and this provides an estimate of the error in the method - we should have 


a value of 1.0 in front of the m term in (41), but instead we get 1.00123, an error of ~ .1%. 


As for the coefficient of the term, it agrees well with the accepted value of —.163 ITS]. 
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B. Self-consistent, Self-Coupled 


We use the same setup and parameters to find the ground state energies for the modified, 
self-sourcing gravity system. In Figure we have the plots corresponding to Figure 
showing the functions P and $ for the smallest and largest mass values. 



FIG. 3. The numerical solutions for P and <1> for m = .7 (top) and m = 2 (bottom) for the 
self-coupled system. Here, i?oo = 200 with AR = 200/5001, and we use e = 10“® as the tolerance 


for the iteration (see (40)). 


This time, we fit to the function: A guided by the form of (24), and obtain 

2.0205 m 

C(m) = 


(42) 

2 +0.45448 m^ ^ ^ 

The numerically-determined ground state values and best fit curve are shown in Figure 
While the fit curve captures the basic behavior of the ground state energies in this case, 
it does not fit as well as in the SN case. The freedom in generating energies using Mp in 
addition to m allows for a more complicated spectrum (for SN, is the only possibility), 
and we do not expect the simple estimate from the Bohr method to work as well. Still, the 
structure of the spectrum is described by that estimate. 


C. Comparison 

We have the numerical spectrum from the self-coupled case, and can compare that with 
the ffi — .163 m® scaling from the SN system - that is shown in Figure [^- for small masses. 
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FIG. 4. Ground state energy as a function of mass for self-coupled gravity - the points 
numerically determined using the numerical method described in the previous section, and 


curve is the best fit curve, C{m), from (42). 


are 

the 


the two ground state energies agree well, but they begin to diverge near rh = 1.2. Note that 
both energies are offset by the same amount (these are bound state energies, so we expect 
them to be negative, they have been shifted upwards by the constant factor mc^, the value 
of the potential energy at inhnity). The energies are off by around 8% by rh = 2. While the 



FIG. 5. Points show the ground state energy for the self-coupled system, while the solid line is the 
SN curve for the ground state energy. 

SN energies become arbitrarily negative as fh increases, the self-coupled spectrum looks like 
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it asymptotically approaches zero (or —mc^, relative to the value of the potential at spatial 
inhnity) - indeed, this is predicted by the Bohr estimate (23), where the energy scales like 


1/m^ for m large. 

Another way to compare the spectrum of the self-coupled case with the SN energies is 
to look at the energy emitted during a transition from the hrst excited state to the ground 
state. In the SN case, the hrst excited state has energy [H |18]: E 2 = m — .0308 
and so the difference between the hrst excited state and the ground state is: 

Asn = .1322 —43 


or AgN = .1322 in our dimensionless variables. 

To compute the hrst excited state using our approach, we simply perform our iteration 
using the eigenvector associated with the second smallest eigenvalue (a simple modihcation 


of (39))- the energy, as a function of m together with the best ht from (24) with n = 2 (so 
that the ht function is Am/{8 + B m^)) is shown in Figure 


E2 



FIG. 6. The energy of the first excited state for the self-coupled case, together with its best-fit 
curve from the Bohr estimate. 

Subtracting the ground state for each mass, we get the Age to compare with the SN case 

- that diherence is shown in Figure together with the Asn from above. Notice there that 
by m = 2, the energy diherence for SN is roughly ten times that of the self-coupled case 

- the frequency of emission (if radiation were emitted) for SN at this mass value would be 
ten times the frequency in the self-coupled case, providing a clear signature for one over the 
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FIG. 7. The energy difference between the first excited state and ground state for SN (solid curve) 
and the self-coupled case (points). 


other. If (23) continues to hold qualitatively for larger n, we see that the self-coupled case 
has a maximum energy of me? (the limit as u —)■ cx)), so any transition energy is bounded 
in this case, while the transition energies of SN will always be of the form: afh^ for some 
constant a. 


D. Total Energy in a Boson Collapse Model 


In the original application of the SN equations, appearing in [T], the total energy of N 
bosons in their ground state was calculated - the Schrodinger-Newton pair appears as 


i h \ 'I'(r,t) 

ot 2 m 

V^<h(r,t) = A ttG N m\b*(r,f) \b(r,f) 


where N particles of mass m are interacting gravitationally in the same state. This is a very 
different application of the SN system (as compared with single-particle collapse), and yet 
we can use our self-consistent scalar gravity in place of the gravitational held equation here, 
just as for the single particle case. We’ll replace the Poisson equation with 


V^<I)(r) = 


47rG 


Nm T*(r, t) ^(r, t) <I>(r) -|- 


2$(r; 


V<I>(r) ■ V$(r). 


as usual. 


(45) 
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The total energy, from the SN approach, is T^tot = Nmc^ — .163 , a dimen¬ 

sionless energy of = rfiN — .163 m® - that is plotted in Figure (the choice of m 
just changes the scale in N, so we have left an unsealed axis there). There are three distinct 
regions in the total energy curve - in the first, where the derivative of the total energy is 
positive, adding particles adds energy. In the second, where the derivative is negative, but 
Etot is still positive, adding particles decreases the total energy. Finally, the total energy 
becomes negative when using the SN equations. 



FIG. 8. The total energy, from [T], for SN applied to a system of N self-gravitating bosons in the 
ground state. 


We can calculate the corresponding total energy using the self-consistent scalar gravity 


in the form (45) ~ taking m = .3 and m = .5, we hnd the total energy and plot them with 
the SN result at those masses in Figure The total energies are different between the two 
cases - the peak that separates the region in which the total energy grows with additional 
particles from the region in which the total energy decreases with additional particles has 
shifted (upward in N for the self-coupled case), and the self-coupled energy does not become 
negative. There is additional physics which must be introduced (and is considered in detail 
in m) as the regimes change, and our goal here is only to once again draw a distinction (in 
principle detectable) between using the Poisson form of Newtonian gravity and using the 
nonlinear, self-coupled scalar gravity. 
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FIG. 9. Total energy for N bosons, each with mass m = .3 (top plot) or m = .5 (bottom plot). 
The points represent the total energy as calculated using self-consistent scalar gravity, while the 
solid curve is the SN prediction from [T]. 


VII. CONCLUSION 

We have introduced an iterative numerical method to find the ground state energies 
for coupled quantum mechanical/gravitational problems. We tested the method with the 
familiar Schrodinger-Newton pair of equations, and found that the method worked well, 
agreeing with known results for that system. In addition, we used the Bohr method to 
predict the form of the ground-state energy for SN as a function of mass. Then, we used the 
same approach for the self-consistent, self-sourced scalar gravity from 1121 . The spectrum 
has a very different mass dependence there, and this is to be expected given the additional 
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(and necessary, for a theory of gravity) nonlinearity appearing in the gravitational field 
equation. Once again, the Bohr method provided a qualitatively relevant estimate for the 
energy dependence on mass for both the ground state and hrst excited state. The difference 
in energies provides a distinct signature for the self-coupled gravity - the transition energy 
from the hrst excited state to the ground state is much less in the self-coupled case than in 
SN. 

It would be interesting to explore the dynamics of the new, self-coupled case, and compare 
with the dynamics of SN, which are known. Aside from specihc computational targets, we 
use special relativity to motivate the self-sourcing in the gravitational held equation, but we 
use the non-relativistic Schrodinger equation to capture the quantum mechanical behavior 
- we could probe the high energy solutions by using the Klein-Gordon or Dirac equations. 
In these cases, it would be interesting to see how the self-coupled gravitational held arises 
in the context of [I9] in which the authors show that Schrodinger-Newton is a natural limit 
of helds coupled to gravity, and the Poisson equation for the gravitational held appears as 
the hrst term in their expansion. 
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